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We study the breathing oscillations in bose-fermi mixtures with Yb isotopes in the largely 
prolate deformed trap, which are realized by Kyoto group. We choose the three combina- 
tions of the Yb isotopes, 170 Yb- 171 Yb , 170 Yb- 173 Yb and 174 Yb- 173 Yb , whose boson-fermion 
interactions are weakly repulsive, strongly attractive and strongly repulsive. The collective 
oscillations in the deformed trap are calculated in the dynamical time- development approach, 
which is formulated with the time-dependent Gross-Pitaevskii and the Vlasov equations. We 
analyze the results in the time-development approach with the intrinsic oscillation modes 
of the deformed system, which are obtained using the scaling method, and show that the 
damping and forced-oscillation effects of the intrinsic modes give time-variation of oscilla- 
tions, especially, in the fermion transverse mode. 



I. INTRODUCTION 



Over the last several years, there have been significant progresses in the production of ultracold 
gases, which realize the Bose-Einstein condensates (BEC) [l- 



, two boson mixtures 

] . In particular the 



degenerate atomic Fermi gases [7j], and Bose- Fermi (BF) mixing gases ABC 
BF mixtures attract physical interest as a typical example in which particles obeying different 
statistics are intermingled. Using the system we have a very big opportunity to get various 
new knowledge about many-body quantum systems because we can make a large variety of 
combinations of atomic species and control the atomic interactions using the Feshbach resonance 



111- 



Theoretical studies of the BF mixtures 
jhase diagram and phase separation 
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26, 



lave been done for static properties 
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171 ]. for the 



1814211]. for stability [22H24I] and for collective excitations 



281434]. 



Recently Kyoto university group has performed researches on the trapped atomic gases of the 
Yb isotopes: the BEC [35J and the Fermi-degeneracy 36]. They also succeeded in realizing 
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the BF mixtures of the isotopes. The Yb consists of many kinds of isotopes — five bosons 
^168, 170, 172, 174, 176 y d ) anc j two fermions ( 171,173 Yb), which give a variety of combinations in the 
BF mixtures. The scattering lengths for the boson-fermion interactions have been determined 
experimentally by the group 



36 38] , and the observation of the ground state properties and the 
collective oscillations of the BF mixtures is now under progressing. 

In many characteristic properties, the spectrum of the collective excitations is an important 
diagnostic signal of these systems; they are common in many-particle systems and are often 
sensitive to the interaction and the structure of the ground and the excited states. 
We have constructed a time-dependent dynamical approach with the time-dependent Gross- 



Pitaevskii (TDGP) and Vlasov equations, and have studied the monopole 



oscillations 



31 



321 ] and dipole 



33] of the BF mixtures. In these works, we have found that collective motions in the 



BF mixtures include various modes such as the boson and fermion intrinsic modes, the forced 
oscillation modes and so on. Then the intrinsic frequencies of these modes in our approach 
are different from those in the sum-rule [3] and the scaling 
approximations of the random phase approximation (RPA) 



39|;42j approaches, which are 
Furthermore the time- 



dependent behaviors in our approach are the same as those in RPA only in the early time, and 
show difference in latter time. As the boson-fermion interaction becomes stronger, the difference 
appears in earlier time. Thus the BF mixing gases show new dynamical properties different from 
the other finite many-body systems such as atomic nuclei. 

We choose three combinations of the Yb isotopes, 170 Yb- m Yb , 170 Yb- 173 Yb and 174 Yb- 173 Yb , 
whose boson-fermion interactions are weakly repulsive, strongly attractive and strongly repul- 
sive. In the previous paper j^ . we have calculated the quadrupole oscillation of these mixtures 
in the spherical trap using the time-dependent dynamical approach, and then compared the 
results with those obtained in RPA. The RPA predicts the same intrinsic modes as the dy- 
namical approach, but it describes the whole oscillations of the mixture as linear combination 
of the intrinsic modes. When the BF interaction is weak, and/or the amplitude is small, the 
linearly-combined oscillations can be good approximation, so that the RPA gives the results 
consistent with those in the dynamical approach. However, when the BF interaction is strong, 
or the amplitude is large, the oscillations of the BF mixtures obtained in the dynamical approach 
show the very different behaviors from the RPA; in the strongly-attractive BF interaction, the 
fermi gas overflows from the boson occupation region and makes large expansion, which leads 
to a monopole oscillation mode. In the case of the strongly-repulsive interaction, the fermion 
intrinsic-mode has rapid damping, so that the fermi gas comes to co-oscillate with the bose gas. 
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In actual experiments, the largely-deformed trap of axial-symmetry is used, and the oscillation 
behaviors may be different from those in the spherical trap. The monopole and quadrupole oscil- 
lations, which are coupled in axial-symmetric traps, are recombined into two modes oscillating 
along the symmetric axial and the transverse directions. Such behaviors are observed in the 
two-component fermi-gas oscillations but no studies have been done for the BF mixture. 
In this paper, we investigate the breathing oscillations of the BF mixtures in the largely deformed 
system of the prolate shape. Then we examine oscillation behaviors along the longitudinal and 
transverse directions separately, which are defined as the directions of the symmetry axis and 
perpendicular to it, respectively. 

In the next section, we explain the dynamical model to calculate the time evolution of the system. 
In Sec. Ill, we examine the intrinsic collective modes and their coupling behaviors in the scaling 
method. In Sec. IV, we show the time-development calculations in the TDHF+Vlasov approach 
for the breathing oscillations in the BF mixture, and discuss their properties. The summary of 
the present work is given in Sec. V. 



II. FORMALISM 
A. Total Hamiltonian 



In this work, we consider the mixture of a dilute boson and one-component-fermion gases at 
zero temperature in the trapping potential of the axial symmetry; the symmetry axis is chosen 
to be the z-axis. We assume zero-range interactions between atoms, and no fermion-fermion 
interactions. The hamiltonian is given by 

h 2 



H 



d 3 q 



+ 2j^asB { ~ Hqm}2 



M B 



r, 2 



2M 



2 ^ a %% +MF) ^q)Hqmq)m 



(1) 



where (ft and tp are boson and fermion fields, M B and Mp are the boson and fermion masses, Qb 
and Vlp are the transverse trapped frequencies of the boson and the fermion, obb and clbf are 
the boson-boson and boson-fermion s-wave scattering lengths. The q = (q?, q£) is the positional 
coordinate, and kl is the trapping-potential frequency ratio of the longitudinal to the transverse 
directions. 
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With the boson harmonic oscillator length for the transverse direction Rb = (H/Mb^b) 1 ^ 2 , we 
introduce the dimensionless coordinate r = (?"i,r2, 7*3) = (rx,z) = q/RB, and he scaled boson 
and fermion fields, <j) = Rb~^^ an< ^ ^ = -^b^V'- Using the variable and fields, the scaled 
hamiltonian H = H/HQb becomes 

H = jd\ -I^t (r) V^( r ) + l( r 2 + 4 z 2^t (r) ^ (r) + £B£ { ^t (r) ^ (r)} 2 

-^^(r)V 2 r i>(r) + -m f uj 2 f (r 2 T + k\z 2 )^ (r)^(r) 



+/iBF0 t ( r )</ , ( T ')V ;t ( r )V ; ( r ) 



(2) 



where rrif and are the fermion-boson ratio of the mass and the trapping-potential frequency, 
rrij = Mf/Mb and ojf = Qp/^B, respectively, and the dimensionless coupling constants of 
the boson-boson and boson-fermion interactions are defined by gsB = ^haBBR^ 1 and /irf = 
AnhaBF'mfO- + m f)~ 1 R'B^ ■ 

At zero-temperature, the total wave function of the BF-mixture including Nf, bosons and Nf 
fermions is written by 



N h 



(3) 



where 4> c is a wave function of the condensed boson and ^/ is a Slater determinant of fermion 
single-particle wave functions ip n . 
The total energy is given by 



Ef 



— E V r ^(r)V r Vn(r) + \m^)(r 2 T + k\z 2 ) ^ ^(rty„(r) 



+h B F(t>l(r)<f> c ( r ) Yl ^n( r )^n(r 



(4) 



B. Time Evolution Equation 



The time evolution of the wave functions are obtained from the variational condition: 



5 j dr^{r)\ (i§^-Hj |$(r)> = 0. 



From the condition, the coupled TDGP and TDHF equations are derived: 



• 9 a ( \ 



lv 2 + U B (r)\ o,.(r.r). 



1 



2m f 



V 2 r + U F (r)\ Vn(r,r). 



(5) 

(6) 
(7) 
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The effective potentials Ub and Up are defined by 

Ub{t) = ^(r?p + k 2 l z 2 ) + g B BPB(r) + h B FPF{r), (8) 
U F {r) = ~m f u){r?t + K 2 L z 2 ) + hBFPB{r), (9) 

where ps and p F are the boson/fermion densities: 

PB {r) = N b \Mr)\ 2 , (10) 

occ 

Pf{t) = J>n(r)| 2 . (11) 

n 

The number of fermion states necessary for calculation are usually too large to solve the above 
TDHF equations directly, so that we use the semi-classical approach. In the semi-classical limit 
(h — » 0), the TDHF equation is equivalent to the Vlasov equation [43j: 

d f d p 



-f(r, p; r) = — + • V r - [V r U F (r)] ■ V p f(r, p; r) = 0, (12) 
dr i<jt rrif J 

where f(r,p;r) is the fermion phase-space distribution function: 

f(r,p,r) = J dM^(r + \u,A^ (r-^u,r\ [<&} e -* u . (13) 
The time-evolution of the system is obtained by solving Eqs. ([6]) and f j 1 2 j> . 

C. Ground State 

In the ground state, the condensed-boson wave function <p c is obtained from the solution of the 
Gross-Pitaevskii (GP) and the Hartree-Fock (HF) equations: 

- l -Vl + U B {r)} Mr,r), = /^ c (r), (14) 

1 -V 2 r + U F (r) \ ip n {r,T) = e n ^ n {r,T), (15) 



2m f 

where pt> is the boson chemical potential, and e n is the fermion single-particle energy of the n-th 
state. 

In the semi-classical limit, the fermion phase-space distribution function in the HF approxima- 
tion becomes the same one in the Thomas-Fermi (TF) approximation: 

f(r,p) = 9[p f -e(r,p)], (16) 

where pf is the fermion chemical potential, and the TF energy e(r,p) is defined by 

e(r,p) = -^-p 2 + U F (r), (17) 
2m f 
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In this TF approximation, the fermion density p F is obtained from the TF equation: 

' {6ir 2 PF (r)} 2/3 + ^m f u}(r 2 T + K 2 L z 2 ) + h BF p B (r) = p f . (18) 



2m/ 2 

When Nb <C 1, in good approximation, the boson density p B is obtained from the TF equation: 

~(r| + k|z 2 ) + g B BPB(r) + h BF p F (r) = p b . (19) 

In actual calculation, we use the boson and fermion densities obtained in the TF equations (|18|) 
and Eq. f)19|) as the initial inputs of the iteration; using the TF densities in Eq. (j!4|) . we obtain 
the boson wave function <p c , which gives the next-step p B in Eq. (jlOp . We iterate solving c/> c and 
searching the fermi energy p,f in Eq. (|18|) to give the correct fermion number, 

III. BREATHING OSCILLATIONS IN THE SCALING METHOD 

The deformed system of BF mixtures generally has the coupled oscillation of the monopole and 
quadrupole modes, which are the breathing oscillation. 

The mode-structure of breathing oscillation in the deformed system has been studied, for exam- 
ple, in two-component fermi gas |42f j . but it is not clear for the deformed BF mixtures. Thus, 
before presenting the calculation of the time-evolution, we should clarify the mode-structure of 
the breathing oscillations in the deformed BF mixture using the scaling method as in Ref. 
it is a macroscopic approximation consistent with the energy- weighted sum rules \^\ - 

A. Scaled Wave Functions 

In the scaling method, the scaled wave functions are introduced on the basis of the ground-state 
one-body wave functions: 

A (r,r) = e^ s(r ' r) e AT - s(r)+ 5 A i>sW^ c(e A T , fl (r) rT . e A i!S (r) z ^ ^ 

where 

Ur, t) = 1 [XtA^t + Ar,«(r)z a ] , (a = B,F). (22) 

The factor exp(z£ a ) is the Gallilei-transformation factor, which is necessary for the scaled wave 
functions to satisfy the continuum equation. 

The variables, \t,b, ^l,b, ^t,f, ^l,f, are the collective coordinates describing the boson longitu- 
dinal breathing (BLB), boson transverse breathing (BTB), fermion longitudinal breathing (FLB) 
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and fermion transverse breathing (FTB) oscillation modes, and A's are the time-derivatives of 
them. In the BF mixtures of g BB = h BF = (no interactions), the matrices B and C become 
diagonal, so that the four breathing modes are completely decoupled. 

Substituting the wave- functions (120|21j) into the total energy functional ([4]), we obtain the total 
energy: 

E T = \f d 3 r{(V r Z B ) 2 p B + e 2XT - B (T Btl +T B:2 ) + e 2XL < B T Bt3 } 

+ 1 J d 3 r{e- 2X ^(r 2 + rl)+e- 2X ^Kiz 2 )}p B (r) 
+ 9j^ e 2X T , B+ x L , B J d 3 r/9 | (r)i 

+ 2^7 / ^ { m K V ^)V + e2XT ' F ( T F,i + T F ,2) + e 2A ^T>, 3 } 
+ ^m f io 2 f J d 3 r{e- 2XF >°(r 2 + r 2 )+e~ 2XF -K 2 L z 2 )}p F (r), 

+ h BF e 2XT ' B+XL > B+2XT > F+XL > F J d 3 r p B (e XT < B r T ;e XL ' B z)p F (e XT ' F r T ;e XL > F z), (23) 



where 

T B 



J d 3 r|V r 0(r)| 2 , (24) 



(z = 1,2,3). (25) 



In order to evaluate the ground state of the system we use the Thomas-Fermi (TF) approx- 
imation, where the momentum distribution becomes spherically-symmetric, the T B and T F i 
become 

T B (r) = 0, (26) 
T>,i(r) = 2>, 2 (r) = T F , 3 (r) = ±(6vr 2 )§ [p F (r)]l . (27) 

In order to simplify the analytic equations, we use the dimensionless variables: 

t 1 h BF m)uj)g BB 

h = 2 ' x= o 2 (n,r2,K L z), (28) 

m fojj g BB Stt z 

2m)ufg BB 2my 2 g 3 BB 

~ ^4 PB, n F ~ ^4 PF> ( 29 ) 

2m 8 f u}°g 2 BB 2m}uj s f g 2 BB 
e B = ^ 1Mb, e F = — 4 up, (30) 

we obtain the scaled total energy, which is also dimensionless: 

, (0 ) 2KLm 2 f 8 uj 3 f 5 g% R 
E t = £Ji BB 3r(\ = 0, A = 0) 

f 1 3 £ ) 

d 3 xlx 2 n B -\ — n 2 B H — rip + x 2 n F + hn B n F > , (31) 
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where x 2 = \x\ 2 . The variations of the total energy, SE^ ' /Sub = and 5E^ /5np = , give 
the TF equations: 

ns + hnp = es — x 2 , (32) 

2 

np + hnp = ep — x 2 . (33) 
From the equation of the second variation: 

we obtain the stability condition of the TF solution: 

UF < 2W' (35) 
B. Collective Oscillations 

Now we study the collective oscillation in the scaling method. The total energy for collective 
oscillations (j3|) is written by 

1 ~ f 2e~ 2 ^ T - s ■ g-2A i S _ 

St = 2 Xb {-^— X Ib + ^-^l,b 

+ J X B + e 2XT ' B+XL ' B V bb 

1~ 2 ( 2e ~ 2XT ' F \2 e~ 2XL < F - 2 

+ 2 AirW / | § T,F + 3k 2 A L,F 

+ ± T F + ± X F + V bf (36) 



where 



fp = ^l d 3 xnp(x) (37) 



X b ,f = j d i xx 2 n B ,F(x) (38) 
Vbb = J d 3 xn 2 B (x) (39) 

2^T,B+^L,B+2^T,F+^L,F 



xfc J A„ B (eW,e^*3)M^ T ,e^* 3 ) (40) 



Using the vector notation for the collective variables, *A = (\t,b, ^l,b, At,f, Al,f)i and expand- 
ing the total energy to the order of 0(A 2 ), we obtain the oscillation energy of the system: 

AE T = E T - EP ph -* ABA + ^*ACA. (41) 
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The matrices B and C are denned by where 



B 



and 



C 



( \X B + W hh 
2V bb 





3 







2V bb 

3^B + V bb 






3k- 



-x B 






3w 






-X, 









x, 










L f 







+ 



I* 

15 



*^3 



T3 



1^1 - 4% 

- S^l - 5^3 



\{f F + X F ) ) 

2_ 

15 



2 v 3 



15 



v 2 - ±v 3 



4% ^ 



where 



Vi 

v 3 



d 3 xxn B 
d 3 x xn F 
d 3 x x 2 



dn F 
dx 
dn B 
dx ' 
driB dn F 



dx dx 



(42) 



(43) 



The stability conditions of the breathing-oscillation states are obtained from dE F /X 
(|4ip. they become 

-2X B + 3V bb -V 1 = 
2T F - 2X F - V 2 = 0. 



(44) 
(45) 
(46) 

a — 0; for 

(47) 
(48) 



Eq. (1411) shows that the classical equation of motion for A is harmonic, so that the oscillation 
frequency uj and the corresponding amplitude A are obtained from the characteristic equation: 



{Buj 2 - C) A = 0. 



(49) 



In the case of no interactions {g BB = h BF = 0), the matrices B and C become diagonal; thus 
the four modes of BLB, BTB, FLB, FTB are completely decoupled. 
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Boson- and Fermion-Intrinsic Oscillations 



The collective oscillation obtained from (|49h describes the mode where bosons and fermions 
oscillate with the same eigen- frequency oj (co-oscillating mode). In RPA or the time-evolution 
approach, the BF mixtures sometimes oscillate in the mode where the bosons and the fermions 
oscillate with the different frequencies. Thus, in the opposite case for co-oscillating modes, we 
should obtain the boson- and fermion-intrinsic oscillation modes where only the bosons or the 
fermions in the BF mixture oscillate collectively, especially in the case of large deformation 
(k l < 1). 

1. boson-intrinsic modes: 

The boson-intrinsic modes are obtained by putting the constraint X a> F = (a = T, L) in (l4*Tj) : 
the characteristic equation become 




(50) 



2X B ( 8 - p B - 2u 2 2 - \v B \ ( X t ,b 

Xl,b 

where vb = V^/Xb- The eigen- frequencies u = ojb are obtained by 

u\ = fl " ^v B ){ 2 + ^4 ± \l 4 - 44 + -k\ \. lol l 




In the spherically-symmetric case (kl = 1), the ojb in (|5ip become uj b = 2(1 — u^/lO) and 
5(1 — ub/10); the corresponding amplitudes are (Xt,b, Xl,b) = (1, —2) (quadrupole oscillation) 
and (\t,b,^L,b) = (1)1) (monopole oscillation) each other. 

In the limit of large deformation (kl "C 1), the boson-intrinsic frequencies ujb become 

u\(l-^v B \ ~4 and \k\. (52) 
For uj b w 4(1 — ws/10), the eigenvector is given by 

^t,b \ / I 

It is just the transverse oscillation (BTB) mode. The eigenvector of the mode with u B 
§4(1-vb/10) is 



Xl,b ) \ 4 




(53) 



A T ,B \ / "I + T?4 \ / "I 



Al r / V 4 + 44 / V 4 



(54) 



17 n L 

5 2 
34 K L 

where the bosons oscillates both the longitudinal and transverse directions in out-of-phase. This 
mode is similar to the quadrupole oscillation, but the transverse-to-longitudinal amplitude ratio 
is different (1:4). According to Ref. [461 ]. we call it the boson axial-breathing (BAB) mode. 
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2. fermion-intrinsic modes: 

Putting X a: B = (a = T, L) in we obtain the characteristic equation for the fermion 

intrinsic mode: 

0, 



2X F 



V 



10 



^3 




(55) 



~5 — T 

K L ] 



where v% = V2/X-F and ^3 = V^/Xp. The eigen- frequencies w = wj? are given by 
2(1 



10J + H + 8 40 j KL 



±2< 



1 



«3 
10 



n , V 2 3v 3 , ,, 
H K 

8 40 



+ I ( '-2 + 



t'3 



(56) 



In the spherically-symmetric case («£ = 1), they become oj' F /u'j = 4 + V2 — and cop/u'j 
4 — U2/2 — U3/2; the corresponding amplitudes are given by (Xt,f, ^l,f) = (lj — 2) (quadrupole 
oscillation) and (\tf,^lf) = (1)1) (monopole oscillation) each other. 
In the limit of kl <C 1, the cj^ becomes 

,2" 



2 / 2 
UJ B /U f 



1 — — ) and 
10/ 



4 1 + 



^2 3^3 \ l(«2 + f 



8 40 



(1-1) 



K 2 

L • 



(57) 



The mode with oj 2 ~ 4(1 — V3/IO) is the transverse oscillation mode: 




1 K 2 




(58) 



We call it the fermion transverse breathing (FTB) mode, for convenience. The eigenvector of 
another mode is given by 




02- 



(1=1) 



(59) 



1 



where the FLB and FTB modes are coupled through the boson- fermion interactions. We call 
this new type of oscillation the fermion axial breathing (FAB) mode. 

In Fig. [TJ we show the frequencies of the monopole (a) and quadrupole (b) oscillations in the 
spherically-symmetric case. The frequencies of the quadrupole oscillation have no level crossing, 
and co-oscillating collective modes (dotted line) are almost the same with the boson-/fermion- 
intrinsic modes. In the monopole oscillation, the frequencies of the boson and fermion modes 
have crossing points at 1ibf/9bb ~ —1.5 and 2. In the region of —1.5 < Iibf/qbe % 2, the 
co-oscillating collective oscillation almost agrees with the intrinsic modes, and, in the region, 
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1ibf/9bb % —1-5 and 2 < Hbf / qbb ■> the co-oscillating mode is found to be different from the 
boson-/fermion- intrinsic modes. 

In the same figure, we plot the boson- and fermion-intrinsic frequencies in the TDGP+Vlasov 
approach (the open squares and the full circles) j^l], Qj- We find that the scaling method well- 
reproduces the results in the TDGP+Vlasov approach for the boson-fermion attractive interac- 
tion, but gives the different results when Hbf / 9bb ^ (monopole oscillation) and Ubf/obb ^ 1 
(quadrupole oscillation). The similar behavior is also found in the case of the dipole oscillation 



331 ] : We have explained on this discrepancy in the previous papers 



0,Q. 



IV. RESULTS AND DISCUSSIONS 



In this section, we show the calculations on the breathing oscillation mode of the deformed BF 
mixtures calculated in the time-development approach. 

Kyoto group controls the interactions by changing the combinations of the Yb-isotopes. We 
search combinations which are available in the present calculation, and choose the three kinds 
of combinations 170 Yb- Yb , 170 Yb- 1 73 Yb and 174 Yb- 173 Yb, where the boson-fermion interac- 
tion is weakly repulsive, strongly attractive and strongly repulsive, respectively. The coupling 
constants are shown in Table HI 

We deal with the Yb-Yb system, where the number of the bosons and the fermions are Nf, = 
10000 and Nf = 1000, respectively. The parameters for the trapping potential is Q,b = 27r x 300 
(Hz), k jj = 1/6 and Uf = 1: the same trap in the actual experiments. Furthermore, we use the 
same mass for all Yb-isotopes (mf = 1); the mass differences can be safely neglected for Yb 
isotopes.. 



A. Ground States 

In Fig. [21 we show the boson and fermion density distributions of the ground states for 170 Yb- 
1 Yb (upper panels), 170 Yb- 173 Yb (middle panels) and 174 Yb- Yb (bottom panels). The 
boson densities are center-peaked in all cases. In contrast, the fermion densities have a center- 
peaked shape largely distributed in the boson-distributed regions in attractive BF interactions 
hsF/gBB < ( 170 Yb- 171 Yband 170 Yb- 173 Yb in Figs. Eh and Ha), and have a surface-peaked 
shape when Iibf/9bb > 1 ( 174 Yb- 173 Ybin Fig. [2b). This BF interaction dependence of the 
fermion density distribution can be easily explained in the TF approximation, where the ground- 
state density is determined not by hsF and gsB but through the ratio Hbf y 1 Qbb- 
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B. Description of Oscillations 



In order to solve the Vlasov equation <\12h numerically, we use the test-particle method |47l |: 
using the test-particle coordinates and momenta and p i: the fermion phase-space distribution 
function is written as 

3 N T N f 

f(r,p,r) = { -S- ^ 6{r - Ti(T)}8{p - pt(r)}. (60) 
where Nt is the number of test-particles per fermion. 

Substituting Eq. (|60p into Eq. (|12p . we can obtain the equations of motion for test-particles: 

4n(r) = ^, (61) 
ctr raj 

± Pi {j) = -V r U F (n). (62) 

As the initial conditions of time-development at r = 0, we use the boosted condensed-boson 
wave function and the fermion test-particle coordinates: 

Mr,r = 0) = exp |l(6 T 4 + 6 L z 2 )| $\r), (63) 

p T (i) = p)f(i) + m f uj f CTr T (i), p z (i) = p[ 9) {i) - 2m f K L uj f c L z(i) (64) 

where br, &l, and are the boost parameters, and the superscript (<?) represents the ground 
state. Then, the initial boson/fermion current densities j b,f{ t ) become 

j B (r,r = 0) = {b T r T + b L zz)p%\r), (65) 
j F (r,T = 0) = uj f (c T r T + c L zz)pf\r). (66) 



Solving Eqs. ([6]), (I6ip and (J62J) with the above initial states, we obtain the time-development of 
the condensed boson wave function and the fermion phase space distribution function. 
In order to examine the aspects of the breathing oscillation, we introduce the quantities: 

Ax L (r; S ) = ^^-l, (67) 



where s = B,F, and Rlpt)(B) and Rl{t) are the root- mean-square radii of the bose and fermi 
density-distributions along longitudinal (transverse) directions. The quantities Ax represent the 
change of distributions in longitudinal and transverse directions. 
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In order to inspect the oscillation modes, we use the strength functions defined by the Fourier 
transform of Axl,t(s) (s = B,F): 



In actual calculation, we take < r < 200 for the integration interval unless otherwise 
noted. The strength functions St,l an d Ct,l exhibit the time-odd and time-even components of 
Axlt(B,F;t), respectively. In the linear-response approximation, only the St,l has strength 
for the present initial conditions: Axl,t = and dAxL^/dr ^ at r = 0. 



First, we present numerical calculations in the TDGP+Vlasov approach with the longitudinally- 
deformed initial condition, bx = ct = and hi, = cl = 0.01, where the boson and fermion gases 
are in-phase. We note that this condition is consistent with the actual experiments by Kyoto 
group. 

In Fig. O we show the time-dependence of Ax L>T for 170 Yb- m Yb, 170 Yb- 173 Yb and 174 Yb- 
173 Yb mixtures. The Axl(B) is found to oscillate monotonously (dotted lines), and show damp- 
ing only in 174 Yb- 173 Yb mixture (Fig. [3k, dotted line) where the BF interaction is strongly 
repulsive. The boson Axt(B) seems to be a superposition of two modes with different periods 
(solid lines); the contribution from the longer-period mode is rather dominant. The BTB mode 
seems to contribute very small except in 170 Yb- 173 Yb . 

The Axl(F) of the 170 Yb- 171 Yb and 174 Yb- Yb mixtures show qualitatively similar behaviors: 
slight beat and damp (dotted lines). In the 170 Yb- 173 Yb mixture, the Axl(F) (Fig. [3H, dotted 
line) is smaller than those in the other two mixtures, and comparable with Axt(F) (Fig. Oi, 
solid line). In this mixture, we see the monotonous increase of the Axt(F), which shows the 
expansion of the fermion gas in time. 

In Figs. S( 170 Yb- 171 Yb and 170 Yb- 173 Yb) and[5] ( 174 Yb- 173 Yb ), we show the strength functions 
for the oscillation. We find that the Sl(B) and St(B) have one large peak at the same frequency 
in all mixtures; The ratios of the peak heights are —Sl(B)/ St(B) ~ 4.4 ( 170 Yb- Yb), 5.3 
( 170 Yb- 173 Yb), 4.2 ( 174 Yb- 173 Yb), which are close to the ratio in the BAB mode obtained in 
the scaling method, -S L {B)/S T {B) = 4 in Eq. (l54j) . 

In the 170 Yb- Yb mixture, the Sl(F) (Figs.H^, solid line) shows one peak at u ~ 0.26, which 
is close to a peak of the Sl(B) (dotted line), while the St{F) has a peak around u) w 1.9 




(69) 



(70) 



C. Longitudinal Breathing Oscillations 
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(Figs. Sla). The 174 Yb- 173 Yb mixture shows the similar behavior except the St{F), which shows 
a large peak at oj ~ 0.26. The existences of large peaks in Sl(F) result in the monotonous 
behaviors of the fermion oscillations in Fig. [3l 

In the 1 74 Yb- 1 73 Yb mixture (Fig. [5]), they show that the BTB and FTB modes exist in the 
oscillation, and that they are even-functions of r. As explained in Ref. [34], these BTB and 
FTB oscillation modes are excited by the fermion-gas expansion seen in Fig. [3l though this 
expansion is very slow and small. 



D. Transverse Breathing Oscillations 

In this subsection, we discuss the oscillations caused by the transversally-deformed initial con- 
ditions, br = ct = 0.1 and bi = cl = (BF-in-phase condition), and br = —ct = 0.1 and 

= cl = (BF-out-of-phase condition). 
In Fig. El we show the time-dependence of Axlt for the BF-in-phase initial condition, and 
we find monotonous oscillation in all mixtures. The amplitude modulation of Axl is seen in 
170 Yb- 173 Yb (solid lines, Fig.[6b,d) in contrast with 170 Yb- 171 Yb and 174 Yb- 173 Yb ; it shows that 
a longitudinal mode with the longer period has been excited and superposed in the 170 Yb- 
1 73 Yb mixture. The Axl,t for the BF-out-of-phase initial condition are shown in Fig. [7J the 
oscillations of Axl,t ar e not monotonous and show beats and dampings. In addition, the 
superposition of the longitudinal modes are seen even in 170 Yb- 171 Yb and 174 Yb- Yb (dotted 
lines), especially, in the boson oscillations. 

In transverse breathing oscillations, we find that the relative phase between the boson and 
fermion oscillations finally become in-phase independently of the initial conditions. It is con- 
firmed in Fig. [U where we plot the Axt in later time 115 < r < 145. We will discuss this 
phenomena in the next subsection. 

In Fig. [91 and 1101 we show the strength functions Sl,t of the oscillations in the 170 Yb- 171 Yb and 
170 Yb- Yb mixtures for the BF-in-phase (left panels) and BF-out-of-phase (right panels) con- 
ditions. The strength functions in 174 Yb- 173 Yb are very similar to those in 170 Yb- 171 Yb , and 
omitted here. For 170 Yb- 173 Yb , we plot both S L (1.1 < u) and C L (w < 0.9) in Fig. [TUJ 
In all mixtures, the St(B) has one sharp peak for both in-phase and out-of-phase conditions; 
in contrast, the St(F) has one sharp peak for the BF-in-phase condition but one sharp-positive 
and broad-negative peaks for the BF-out-of-phase conditions (panels c,d in Figs. [9l and 1 1Q[) . The 
broad peaks for the BF out-of-phase condition correspond to the intrinsic FTB modes, but, for 
the in-phase initial conditions, the intrinsic FTB modes do not appear; the boson-gas oscillations 
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seem to influence on the fermion transverse breathing oscillation. 

In the 170 Yb- 173 Yb mixture, the Cl(B) has one sharp peak for both BF-in-phase and BF-out- 
of-phase conditions, and the Sl have very small peaks at the frequency of the peak position 
of St(B). The long-period time-even modes in Cl(B) are considered to be excited through 
the intrinsic BAB modes because they include both the longitudinal and transverse oscillation 
components (|54p . Thus, we find that the longitudinal and transversal oscillations are approxi- 
mately time-even and time-odd in 170 Yb- 173 Yb . It explains that the longitudinal oscillation is 
not excited by the initial condition, and that the oscillation in 170 Yb- 173 Yb cannot be described 
with the superposition of the intrinsic oscillation. 

In order to examine the above-mentioned results, we perform the calculation with the initial 
condition: bp = 0.1 and ct = bp = cl = 0, where only the boson gas is boosted along the 
transverse direction (Fig. [TT]) . In the first one period, the relative phase between the boson and 
fermion oscillations are in-phase in 170 Yb- 173 Yb (repulsive BF interaction) and out-of-phase in 
170 Yb- 171 Yb and 174 Yb- 173 Yb (attractive BF interaction); it is consistent with the spherical-trap 



calculations 



33, 



34f ] . The period of the boson oscillation is almost constant while that of the 



fermion oscillation varies in time; finally, the periods of them come to be equal and the relative 
phase of the boson and fermion oscillations become in-phase. This behavior is the same one 
with that in Fig. [8j 

Now, we give a short explanation for the convergence to the BF in-phase oscillation in the 
transverse breathing oscillations independently of the initial condition (Figs. [8j and [TT]) . As 



shown in Sec lII CI the ground-state boson density is a function of yVj, + k 2 l z 2 in the TF 
approximation; then, the longitudinal component of the effective fermion potential Up in Eq. 
becomes z ■ V t Uf ~ k-l^t • V r t/jr. It shows that, when kl <C 1, the force for the fermions along 
the longitudinal direction should be weak in comparison with that along the transverse direction; 
the intrinsic longitudinal motion of the fermion gas is very slow, and should be damped easily 
in early stage of oscillation. Then, it remains the forced oscillations by the boson gas with the 
longer period, and relative phase between the boson and fermion gases becomes in-phase. 



E. Comparison with the Scaling Method 

In this subsection, we discuss the BF-interaction dependence of the intrinsic-mode frequencies 
obtained in the TDGP+Vlasov approach, and compare them with those in the scaling method. 
Fig. [121 shows the /isi?/g , Bs-dependence of the boson and fermion intrinsic frequencies, which 
are obtained from the peak frequencies of the strength functions. The boson and longitudinal 
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fermion oscillations are found to have small /is^/g^-dependence and the scaling method well 
reproduces the time-dependent calculations for these oscillations. 

The fermion transverse oscillation has large interaction dependence, and the scaling method 
cannot reproduce the TDGP+Vlasov calculation (Fig. [T2l) . especially at h^F / 9bb = 1-325 
(large repulsive BF interaction). In the scaling method calculation, the frequency of the FTB 
mode has the minimum around Hef/qbs ~ 1; and increases with increase of Jt,bf / 9bb m the 
region of Hbf / 1 9BB > 1 (Fig. fT2l solid line); it comes from the structure change of the fermion 
ground-state density distribution into the surface-peaked shape (Fig. [2]). On the other hand, 
in the TDGP+Vlasov approach, the density distribution and the velocity field varies in time 
in large amplitude oscillations. In the mixture of Hbf/9bb > 1> the minimum position of the 
fermion potential is on the spherical surface at the border of the boson density distribution, 
and the fermion gas oscillates around this surface; it results in the frequency of the fermion 

he same discrepancy appears in the dipole 
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transverse mode smaller than that in the sum-ru 
and quadrupole oscillations in the spherical trap 
The longitudinal and transverse modes are found to decouple, and the frequencies of the boson 
and fermion intrinsic oscillations are very close. Thus the boson and fermion gases can easily 
oscillate with the same period; it also explains the convergence into the BF in-phase mode in 
the transverse oscillation, which has been discussed in the previous subsection. 



V. SUMMARY 



In this paper, we investigated the collective breathing oscillation of the BF mixtures of Yb 
isotopes in the deformed trap (kl = 1/6). Three combinations were examined: 170 Yb- 171 Yb, 
170 Yb- 173 Yb and 174 Yb- 173 Yb , where the BF interactions are weakly repulsive, strongly attrac- 
tive and strongly repulsive. 

In the case of the weak BF interactions, the longitudinal and transversal breathing oscillations 
decouple in the largely deformed trap of prolate shape. Because of the large interaction energy 
coming from the condensed bosons, the boson breathing oscillations couple into the intrinsic 
BTB and BAB modes, and the intrinsic FTB and FAB modes are made from the fermion 
oscillations. The actual time-development processes are described by the superposition of these 
intrinsic modes. 

The intrinsic transverse modes (BTB and FTB) have very close frequencies, and it is similar in 
the longitudinal modes (BAB &: FLB); The frequencies of the two pairs separate very well. 
In the spherical trap, the fermion oscillation have been shown to have the boson-forced oscillation 
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modes and the two intrinsic modes which correspond to the inside- and outside-fermion oscilla- 
tions for the boson-distributed regions. In the present prolate deformed trap, no outside-fermion 
modes appear. 

The longitudinal oscillations are rather simple; they are approximately monotonous in the boson 
and fermion oscillations in the all BF mixtures. For the in-phase initial condition, the BAB mode 
and the fermion intrinsic modes appear in the boson and fermion oscillations, respectively. 
The transverse breathing oscillations have complex structures as clarified in IV D, and they 
depend on the strength of the BF interactions. In the fermion transverse breathing oscillations, 
we found that the fermion-gas frequencies vary in time and converges into the same frequencies 
with the boson oscillations and the relative phases of the boson and fermion oscillations become 
in-phase. This phenomenon are shown to be explained by the damping of the intrinsic fermion 
modes and the boson-forced oscillations. 

In this paper, we discussed the collective oscillations of the BF mixtures of Yb isotopes at T = 0. 
In the actual experiments, which should be done at very low but T^O temperatures, thermal 
bosons should give some contributions: the introduction of such effects into the present approach 



through two-body collision terms 



481, 



49j should be done also in the future. 
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System 


9BB 


hBF 


h]3F / 9BB 


170 Yb _171 Yb 


5.598 xl0~ 2 


3.207 xl0~ 2 


0.573 


170 Yb _173 Yb 


5.598 xl0~ 2 


-7.126 xl0~ 2 


-1.273 


174 Yb _173 Yb 


8.902 xlO" 2 


11.799 xl0~ 2 


1.325 



TABLE I: The coupling constants gsB and hsF of the boson-boson and the boson- fermion interactions 
in 170 Yb- m Yb, 170 Yb- m Yband 174 Yb- 173 Yb . 
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FIG. 1: (Color online) The frequencies of the monopole (a) and the quadrupole oscillations (b) in the 
spherically-symmetric system (/tj, = 1) in the scaling method: the co-oscillating (dotted lines), the boson- 
intrinsic (dotted lines), and the fermion- intrinsic (solid lines) modes. The open squares and full circles 
denote the numerical results obtained in the TDGP+Vlasov calculation for the boson- and fermion- 
intrinsic frequencies, which are taken from Refs. 31 j (monopole) and [sj] (quadrupole). 
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FIG. 2: (Color online) The ground-state density distributions of the BF mixtures 170 Yb— 171 Yb (upper 
panels), 170 Yb— 173 Yb (middle panels) and 174 Yb— 173 Yb (bottom panels) at z — (transverse direction, 
left panels) and at tt = (longitudinal direction, right panels) . Dotted and solid lines are for the bosons 
and fermions, respectively. The thin long dotted lines indicate the boson density in the TF approximation. 
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FIG. 3: (Color online) Time evolution of Axl,t{B) (boson oscillation, left figures) and Axl,t(F) (fermion 
oscillation, right figures) for the in-phase initial condition: bx — ct — and — = 0.01. The upper, 
middle, and lower panels are for the 170 Yb- 171 Yb, 170 Yb- 173 Yb and 174 Yb- 173 Yb mixtures. The dotted 
and solid lines represent Axl and Axt, respectively. 
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FIG. 4: (Color online) The strength functions Sl (longitudinal, upper figures) and St (transversal, lower 
figures) of the boson and fermion oscillations for the in-phase initial condition. The left and right panels 
are for the 170 Yb- 171 Yb and 174 Yb- 173 Yb mixtures. 
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FIG. 5: (Color online) The strength functions Sl (longitudinal, upper figure) and St & Ct (transversal, 
lower figure) of the boson (dotted lines) and fermion (solid lines) oscillations in the 170 Yb- 173 Yb mixture 
for the in-phase initial condition. In the lower panel, St and Ct are plotted in u) < 1.3 and 1.4 < lu, 
respectively. 
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FIG. 6: (Color online) Time evolution of Axl,t{B) (boson oscillation, left figures) and Axl,t(F) (fermion 
oscillation, right figures) for the in-phase initial condition: bx = cy = 0.1 and = cl = 0. The upper, 
middle, and lower panels arc for the 170 Yb- 171 Yb, 170 Yb- 173 Yb and 174 Yb- 173 Yb mixtures. Dotted 
and solid lines represent Ax^ and Axt, respectively. 
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FIG. 7: (Color online) Time evolution of Axl,t{B) (boson oscillation, left figures) and Axl,t(F) (fermion 
oscillation, right figures) for the out-of-phase initial condition: bx — ~ct — 0.1 and 6l = cl = 0. The 
upper, middle, and lower panels are for the 170 Yb— 171 Yb, 170 Yb— 173 Yb and 174 Yb— 173 Yb mixtures. 
Dotted and solid lines represent Axl and Axt, respectively. 
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FIG. 8: (Color online) The time evolutions of Axt in later time 115 < t < 145 , for the BF-in- 
phase (left panels) and out-of-phase (right panels) initial conditions in the 170 Yb— 171 Yb (upper panels), 
170 Yb— 171 Yb (middle panels) and 170 Yb— 171 Yb (lower panels) mixtures. Dotted and solid lines represent 
the results of the boson and fermion oscillations. 
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FIG. 9: (Color online) The strength functions Sl (longitudinal, upper figures) and St (transversal, 
lower figures) of the boson and fermion oscillations in 170 Yb— 171 Yb for the transversally-deformed initial 
condition: the left and right panels are for the BF-in-phase and BF-out-of-phase conditions. Dotted and 
solid lines represent the results of the boson and fermion oscillations, respectively. 
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FIG. 10: (Color online) The strength functions Sl (longitudinal, upper figures) and St (transversal, 
lower figures) of the boson and fermion oscillations in 170 Yb— 173 Yb for the transversally-deformed initial 
condition: the left and right panels are for the BF-in-phase and BF-out-of-phase conditions. Dotted and 
solid lines represent the results of the boson and fermion oscillations, respectively. 
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FIG. 11: (Color online) Time evolution of Axt in the 170 Yb- m Yb (a), 170 Yb- 173 Yb (b) and 174 Yb- 
173 Yb (c) mixtures for the initial condition: bx = 0.1 and ct — &l = cl = 0.01. Solid and dotted lines 
are for the bosons and fcrmions, respectively. 
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FIG. 12: (Color online) Intrinsic frequencies of the breathing oscillations of the boson and fermion gases 
obtained in the TDGP+Vlasov approach; circles and squares are for the longitudinal and transverse 
oscillation modes in the 170 Yb— 171 Yb, 170 Yb— 173 Yb and 174 Yb— 173 Yb mixtures Dotted and solid lines 
represent the results of the longitudinal and transverse modes in the scaling method. 



